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ABSTRACT 


The present investigation deals with the postbuckling 
response of the simply supported angle- and cross-ply laminated 
plates under thermal loads. The Rayleigh Ritz energy technique is 
used to obtain the governing nonlinear algebraic equations by 
assuming a double series form of trial functions for the 
displacement fields. A computer program is developed to solve 
nonlinear algebraic equations. The effect of boundary conditions, 
the fiber orientation, the number of lay-ups and the aspect ratio 
on the transverse deflection is studied The analysis is carried 
out for uniform temperature throughout the plate as well as for 
slight temperature gradient through- the-thickness of the plate. 
Results are presented in the graphical form for antisymmetric 
angle-ply and symmetric and antisymmetric cross-ply laminates. 



CHAPTER 1 


INTRODUCTION 

1.1 INTRODUCTORY REMARKS : 

It is well known that the measured strengths of most of the 
materials are found to be much smaller than their theoretical 
strengths. The discrepancy in strength values is believed to be 
due to the presence of imperfections or inherent flaws in the 
material. The flaws in the form of cracks that lie perpendicular 
to the direction of applied loads are particularly detrimental to 
the strength. Therefore, compared with the strength of bulk 
material, man-made filaments or fibers exhibit much higher 
strength along their lengths since large flaws, which may be 
present in the bulk material, are minimized because of small 
cross-sectional area of the fibers. Fibers, because of their 
small cross-sectional dimensions, are not directly usable in the 
engineering applications. They are, therefore, embedded in the 
matrix materials to form fibrous composites. The matrix serves to 
bind the fibers together, transfer loads to the fibers and protect 
them against environmental attacks and damages due to handling. 

Laminated composite plates are made of plies, each ply being 
composed of straight, parallel fibers (e.g., glass, boron, 
graphite) embedded in and bonded together by a matrix material 
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(e.g epoxy resin) Each ply may be considered as homogeneous, 
orthotropic material having a value of Young’s modulus(E) 
considerably greater in the longitudinal direction than in the 
transverse direction. Adjacent plies will have longitudinal axes 
usually not parallel. Cross-ply laminated plates arise in the 
special case when the longitudinal axes of the adjacent plies are 
perpendicular, whereas angle-ply laminates occur when adjacent 
layers are alternatively oriented at angles of +0 and -0 with 
respect to the edges of the plate. 

Fiber-reinforced composite laminates have high 
strength-to-weight and stiffness-to-weight ratios and are becoming 
important in weight-sensitive applications such as spacecraft and 
space vehicles. A re-entry space vehicle encounters hyperthermal 
loading conditions. In certain cases, the thermal load turns out 
to be the primary one and controls the design. The availability 
of composites and use of optimization techniques in the design for 
minimum weight results in structures which are very thin. The 
plate elements, which are extensively used in aircraft 
construction, are subjected to compressive stresses due to 
boundary constraints. Therefore, the thermal loads may cause 
buckling, especially in thin-walled structures. Though the 
investigation of thermally-induced buckling of composite plate is 
in itself an important aspect, the study of the postbuckling 
response, specially in the presence of a temperature gradient 
through- the-thickness of plate, assumes a far greater importance. 
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As with the mechanical loading, composite plate should be expected 
to operate in the postbuckling range with thermal loads as well, 
if structural efficiencies are to be high Restricting the 
temperature to remain below the buckling temperature, or designing 
stiffnesses so that plates do not buckle, are inefficient 
approaches to thermal structural design. Moreover, buckling is 
actually an idealization of the problem. Generally, the geometric 
imperfections in the plate, the support conditions and even 
imperfections in the assumed temperature field preclude the 
occurrence of the ideal buckling (bifurcation), thus necessitating 
the study of the postbuckling response of such plates This forms 
the scope of the present investigation. 

1.2 LITERATURE REVIEW : 

A recent review article by Tauchert (1991) on 
thermally- induced flexure, buckling and vibration of plates gives 
an up-to-date account of buckling of composite plates, under 
thermal loads. Postbuckling of plates is an important topic and 
has attracted the attention of many researches over the past 
several decades. A detailed literature review on thermal 
postbuckling of isotropic and orthotropic plates has been covered 
in the article by Tauchert (1991). 

Most of the early investigations on the subject of thermal 
postbuckling have been devoted to isotropic thin plates and 
shells. Based on the large deflection theory, Gossard (1952) used 
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the Rayleigh-Ritz method to calculate the deflections of flat or 
initially imperfect plates subjected to nonuniform temperature 
distributions The postbuckling behaviour of rectangular plates 
which are heated symmetrically about the two center lines was 
investigated by Forray and Newman (1962a). Soon afterwards they 
solved the same problem by the combined Rayleigh-Ritz and Galerkin 
procedure which permits more flexibility in the choice of boundary 
conditions (1962b). A method for the analysis of a rectangular 
plate heated nonunlformly and allowed to undergo large deflections 
was presented by Wilcox and Clemmer (1964). Mahayani (1966) 
studied the thermal buckling and postbuckling response of a simply 
supported cylindrical panels subjected to parabolic temperature 
distribution along its radial direction. Thermal postbuckling 
characteristics of skew plates with restrained edges were analyzed 
by Prabhu and Druvasula (1976). Following Berger’s method, the 
large deflection of a heated cylindrical shell was investigated by 
Biswas (1978a). He also presented an analysis for the large 
deflections of a heated orthotropic plate (1978b). Banerjee and 
Datta (1979) studied the large deflections of simply supported 
isotropic plates of any given shape under a varying temperature 
distribution. Biswas (1981) obtained an approximate solution of 
the governing equations for a heated orthotropic plate. A simple 
thermal postbuckling analysis of isotropic circular plates was 
carried out by Raju and Rao (1984). 


Compared to isotropic cases, relatively few investigations 
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have been made on the thermal postbuckling behaviour of laminated 
composite plates. Thermal buckling and postbuckling behaviour of 
antisymmetric angle-ply plates, resting on simple supports and 
subjected to uniform temperature rise, was examined by Tauchert 
and Huang (1986). Since the closed-form solutions to the 
resulting coupled nonlinear equations are not easy, a direct 
application of the principle of minimum potential energy had been 
applied to obtain approximate solutions Numerical results for 
graphite/epoxy laminates illustrated the effects of number of 
lay-ups, the ply angle and the aspect ratio of the plate upon 
thermoelastic response Huang and Tauchert (1988a) again analyzed 
the same problem, but here the emphasis was given to secondary 
buckling. The effect of fiber orientation, number of lay-ups and 
the aspect ratio was also carried out. The scope of this study 
was enlarged in a subsequent investigation (Huang and Tauchert, 
1988b) in which the nonuniform temperature field was considered. 
They showed that unlike the uniform temperature problem in which 
the simply supported laminate remains flat until bifurcation 
occurs, under a nonhomogeneous temperature field, the transverse 
deflections occur immediately upon loading In such cases a 
bifurcation may occur, however, at a later stage of thermal 
loading. A major effort in dealing with the problem of 
thermally- induced large deflections by finite element method was 
made by Chen and Chen (1989). They carried out the study of 
geometrically nonlinear response of composite plates under a 
nonuniform temperature field. The analysis showed that the 
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postbuckling behaviour is influenced by the lamination angle, the 
plate aspect ratio, the modulus ratio and the number of layers. 
The large-deflection response of laminates under thermal loading, 
including flat plates and cylindrical panels, was investigated by 
Huang and Tauchert (1990). The finite element formulation based 
on the first order shear deformation theory was employed to 
predict the shell response. The strength limit of laminates had 
been predicted using Tsai-Wu criterion. With change in 
temperature, thermal and elastic properties (e.g. coefficient of 
thermal expansion, El, Et, etc) do not remain constant. Chen and 
Chen (1991) examined the effect of temperature dependent 
properties on thermal postbuckling behaviour of laminated plates. 
It was found that the postbuckling strength reduced significantly 
when temperature-dependent properties were taken into 
consideration. Meyer and Heyer (1992) analyzed the 

thermally-induced geometrically nonlinear response of 
symmetrically laminated plates. The plate response was considered 
due to a temperature increase that was uniform in the plane of 
plate but had a slight temperature gradient through-the- thickness 
of plate. The case of a completely uniform temperature increase 
but with initial out-of-plane imperfection was also considered. 
The study showed that difference between buckling temperatures for 
the fixed and the sliding simple supports is small and that the 
support conditions can have large influence on the 

thermally-induced nonlinear response. It was also shown that the 
geometric as well as temperature Imperfections play an important 
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role in the initial out-of-plane deflections 

1.3 AIM OF INVESTIGATION . 

The forgoing literature review clearly shows that only a few 
research works have been carried out on laminated composite plates 
subjected to temperature loading. It is also clear that no 
results are available for thermally induced large deflections of 
cross-ply laminated plates. The aim of the present work is to 
study the postbuckling response of antisymmetric angle-ply and 
cross-ply laminated plates, under two types of temperature 
fields : (i) uniform temperature increase throughout the plate, 

(ii) temperature increase that is uniform in the plane of plate 
but has a slight gradient through-the-thickness. 

The effect of the number of lay-ups, the fiber orientation 
and the aspect ratio of the plate on thermally- induced nonlinear 
response is considered. 

1 4 STATEMENT OF THE PROBLEM . 


A plate as shown in the Fig. 1.1-1 2, with length a, width b 
and thickness H is considered The reference co-ordinate system 
has its origin at one corner of the plate. The x-y plane 
coincides with the geometric midplane of the plate, and z— axis 
points downwards. The geometric midplane is the reference surface 
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of the plate. The temperature change, AT, is given by 

AT = T(1.0 + c.z/H) (1.1) 

where T is the temperature change at the midplane of the plate. 
The parameter c controls the magnitude of the 
through- the-thickness gradient. Referring to the co-ordinate 
system mentioned, the temperature change at the top of plate is 
T(1.0 - c/2.0), and that at the bottom is T(1.0 + c/2.0). This 
can as well represent an imperfection in the temperature 
distribution which will cause the plate to begin to deflect out of 
plane as soon as the temperature increases relative to the 
reference temperature. 

The problem can be also analyzed for following cases : 

(i) c = 0, i.e , the plate is subjected to uniform temperature 

(ii) Linearly varying transverse temperature distribution given by 

AT « (T.c).z/H where, now, T.c denotes the temperature difference 
across the plate thickness. Cases (i) and (ii) can in fact be 
superimposed to get a distribution of the type given by Eq. (1.1). 
Both the angle-ply and the cross-ply laminated plates are included 
in the study. Two boundary conditions have been considered : 
(i) fixed simple support conditions : all four edges are fixed 
against in-plane normal and tangential displacements as well as 
out-of -plane displacements (Fig. 1.1), (ii) sliding simple supports 
: all four edges are fixed against in-plane normal and 

out-of-plane displacements, but are free to move or slide 
tangentially(Fig. 1.2) . Sliding is assumed frictionless. These 
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boundary conditions can be expressed as : 


Fixed simple supports (referred to as type 1) 


Along x = 0, a 

uo = 0, vo = 0, w = 0, M =0 

X 

Along y = 0, b 

uo ■= 0, vo = 0, w = 0, M = 0 

y 

Sliding simple supports (referred to as type 2) 

Along x = 0, a 

uo = 0, N = 0, w = 0, M =0 
xy i x 

Along y = 0, b 

N = 0, Vo = 0, w = 0, M = 0 
xy y 


( 1 . 2 ) 


(1.3) 







CHAPTER 2 


FORMULATION AND SOLUTION TECHNIQUE 

2.1 INTRODUCTION : 

In most practical applications, the thickness of a plate is 
small in comparison with its smallest lateral dimension, and hence 
Kirchhoff’s hypothesis may be assumed to be valid. This implies 
that tractions on the surfaces parallel to reference plane are 
negligibly small as compared with the inplane stresses, and that 
the inplane displacements are linear functions of z. 

The determination of the buckling temperature load for a 
particular structure, say a thin plate, is essentially an 
eigen-value problem based on linear ( small deflection theory ) 
but this type of analysis yields only the buckling temperature and 
does not give any information about the magnitude of the lateral 
deflection at the onset of buckling (of plate) and beyond. 
Therefore, the use of non-linear equation, based on the large 
deflection theory, becomes essential for studying the behaviour of 
thin plate in the postbuckling regime. The von karman equations 
are most widely used equations for this purpose. Although these 
equations are available in various forms their analytical solution 
is not possible in most of the cases. Out of the many numerical 
methods available to solve these nonlinear equations the Rayleigh 
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Ritz method turns out to be quite convenient for a wide class of 
problems. 


The Rayleigh Ritz technique or simply the Ritz technique is a 

procedure for applying the principle of minimum potential energy 

to obtain approximate solutions of elastic problems. In case of 

plate problems where the three displacement components u, v and w 

are to be determined, the procedure consists essentially in 

assuming that the desired extremal of a given problem can be 

approximated by linear combinations of suitably chosen functions 

in the following form 

x J 

u = V W u (x,y) 

^ u ij ij ,J 
1=1 j=i 


x J 

v = £ E Bv (x,y) (2 D 

1=1 j=i J J 


M 

W = E 
m= 1 


N 

EC w (x.y) 
mn mn 

n=l 


In these expressions A , £ , and C are the undetermined 

1 J i J nrn 

parameters and u ,v and w are the known functions, each of 
r 1 J 1 J mn 

which can be written as F(x) G(y). These functions are chosen 

such that they are continuous, satisfy at least the prescribed 

geometric or kinematic boundary conditions for any choice of these 

parameters, and are capable of representing the deformed shape of 

the plate problem. When Eq.(2.1) is substituted into the 

potential energy expression, the latter becomes the function of 

parameters A , B and C . These parameters are determined from 

r i j i j mn 
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the condition that the total potential energy of the system is 
minimum with respect to them. This condition gives the same 
number of equations for A , B and C as the number of 

1 j i J mn 

parameters taken. The accuracy of the method depends upon the 

choice of the number of parameters and the approximating shape 

functions u ,v and w 

1 J i J mn 

2 2 DISPLACEMENTS AND STRAINS : 

In view of the Kirchhoff ’ s assumption of thin plate theory, 
the inplane displacements u, v and w at any arbitrary point of the 
plate in the x-, y- and z-directions can be approximated by 

u(x,y,z) = uo(x,y) - z.w.x 

v(x,y,z) = vo(x,y) - z.w.y (2.2) 

w(x, y, z) = w(x,y) 

where uo, vo are the mid-surface inplane displacements and a comma 
signifies the derivative with respect to suffix that follows. 

In view of the von Karman theory of large deflection of 
plates, the strain displacement relationship can be written as 

e = 3uo/3x +( 1 / 2 ) (3w/3x) 2 - z (3w/3x 2 ) 

X 

e = 3vo/3y +( 1 / 2 ) (3w/3y) 2 - z (3w/3y 2 ) 
y 

y = Svo/3x + 3uo/3y + (3w/3x) (3w/3y) - 2 z (3w/3x3y) 
xy 


(2.3) 



14 


In short, these can be written as 

[e] = [e°] + [e° L ] + z [K] (2 4) 

where [c 0 ] are the mid-plane strains, which include terms which 
account for stretching due to transverse displacement w. These 
terms are seen to be nonlinear [K] are the curvature changes of 
the mid-plane during deformations These are functions of x and y 
only. 

2 3 TOTAL POTENTIAL ENERGY : 

The total potential energy V of a structure is the sum of the 
strain energy U stored in the rectangular plate and the potential 
energy due to the thernu l forces and moments, i e. 

V = U + V (2.5) 

T 


The strain energy U stored in a rectangular laminated plate, which 
is assumed to be in a state of plane stress, can be expressed as 



( 2 . 6 ) 


Using the Eqs. (A1 . 3) , (A1 . 6) and (A1.10) from Appendix 1 and the 
subsequent integration witl respect to z yields 



b a 

u = 1/2 I I [ A 11 ( u * x + 1/2 w ^) 2 + 2 A 12 ( v >y + 1/2 wfy)* 

o o 
2 

( U,x + 1/2 W,x) + 2 A j6 (u,y + V , x + W,x W , y ) ( U,x + 1/2 wfx) + 

2 A 26 (U,y ♦ V,x + W,x W,y)( V.y + 1/2 wfy) + A ^ V,y + 1/2 wfy) 2 

+ A fi6 (U,y + V , x + W,x W,y) 2 - 2 { B n ( U, x + 1/2 wfx) W.xx + 

2 

B 12 ^ V,y + 1/2 w *y) W.xx + B lfi (u,y + V , x + W, X W, y ) W.xx + 

2 

B 12 ^ U,X + 1/2 W,X ) W,yy + B 22 ( V,y + i/2 wfy) W,yy + 

B (U,y + V,x + W , x W,y)w,yy + 2 B (u, x + 1/2 wfx)w,xy + 

2 B^Cv.y + 1/2 wfy)w,xy + 

2B 66 ^ U,y + V,X + W,X W > y ) W,xy ) + D la w f xx + 2 D W,yy W.xx + 

4 D 16 W.xy W.xx + D 22 W , y y + 4 W.xy w,yy + D wfxy ] dx dy 

(2 7) 

In the preceding expression the subscript o of uo and vo have been 
dropped for convenience. 

The potential energy due to thermal forces and moments can be 
written as 

b a 

V T = J J[-N*( u.x + 1/2 wfx ) ~ N* ( v.y + 1/2 wfy) - 

T 0 0 T T T 

N (u,y + V,x + W,x W,y) + M W,xx + M W,yy + 2 M w,xy ] dx dy 
xy x y xy" 

(2.8) 

The relevant expressions for thermal forces and moments are 
presented in the Appendix 1 (see A. 14). 
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2.4 APPLICATION OF RAYLEIGH RITZ METHOD 


The geometric boundary conditions of the present problem can 
be satisfied by assuming a double series form of a trial function 
for the displacements. They are given below for the two types of 
edge conditions considered (section 1.4). 


Type 1 (i.e. fixed simple support ) 

i J 

uo (x,y) = Y Y u sin(i7ix/a) sin(jTty/b) 
i J 

i=l J=1 
I J 

vo (x,y) = £ E v . . sindmx/a) sln(jrty/b) (2.9) 

1=1 j=i lJ 

M N 

w (x,y) = E Ew sin(mnx/a) sin(n7iy/b) 

mn 

mwl n= 1 


Type 2 (i.e. sliding simple support) 


i j 

uo (x,y) = T V u , , sln(inx/a) cos(jny/b) 

ij 

i =i j=o 


i j 

vo (x,y) » £ £ v cosdmx/a) sin(jny/b) 

l=o j = i iJ 


( 2 . 10 ) 


M N 

w (x,y) = Y F w sin(mirx/a) sin(mry/b) 
mn 

m= 1 n= 1 


These series are substituted in the total potential energy V 
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obtained by using Eqs (2.7) and (2 8) in Eq (2.5) Then the 
Ritz equations are obtained by setting 


av/au^ = o 

av/av^ = o (2.1D 

av/aw = o 

mn 

This results in a coupled set of nonlinear algebraic equations. 
Which can be written in following general form : 


p 0 i J p 0 i 

znz u fi+rzz 
p q i J l J p 


H 


P Q M N 

V f2 + Z Z Z Z W f3 = -g„ 
ij pqmnmn 1 


+ N1 


P 0 I J 

£ £ £ 
p 


H u 


i j 


f 4 


P Q X J 


P Q M N 


+ ZZZ?V ,f5 + Z Z Z Z w 
P q i 


J U 


p q m n 


f 6 = -g + N2 


P 0 I 

ZEE 
p q 


H u 


i J 


(f7 * s 3 > * £ I H V U (fa 


P Q M N 

+ g)+E£2£w f 9 
4 p q m n mn 




+ N3 

(2.12) 


where Nl, N2 and N3 are 
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ba 

Nl=[J(N*(pn/a) cos(pnx/a) sinCqny/b) + 

O O 

T 

N xy qir ' sin(p7ix/a) cos (qrcy/b) ldxdy 

ba 

N2=JJ[N^(q7i/b) sin(prcx/a) cos(q 7 iy/b) + 

oo 

N 3iP n/a) cos(pnx/a) sintqny/b) ]dxdy 

x y 

ba 

N3=J|[M^(prt/a) 2 sin(p7rx/a) sinCqny/b) + 

O O 

M^(qit/b) 2 sin(p7tx/a) sin(q7ry/b) 

2M T (pn/a) (qre/b) cos(p7rx/a) costqrr/b) ]dx dy 

x y 


In Eq. (2.12) fi to f9 are linear expressions, while g^o g g are 

nonlinear expressions. Where g contains cubic as well as square 

s 

terms. To solve these coupled set of nonlinear algebraic 
equations, a computer program has been developed. This program 
uses standard subroutines available in the NAG library for single 
integration, double integration and for the solution of nonlinear 
equations. The solution of various cases are presented and 
discussed in the next chapter . 


i 



CHAPTER 3 


NUMERICAL RESULTS AND DISCUSSION 

3 1 INTRODUCTION : 

The formulation developed In the chapter 2 has been first 
checked by comparing the results with some of those available in 
the literature for antisymmetric angle-ply laminates and 
symmetric laminates. It is then employed to estimate the 
transverse deflection of antisymmetric angle-ply and symmetric 
and antisymmetric cross-ply laminates, each consisting of N 
equal -thickness layers. The a/H ratio (plate length to 
thickness ratio) considered Is 150. The material properties 
used in the analysis represent a graphite-reinforced composite. 
These properties are : 

E « 155 GPA E - 8.07 GPA G lt = 4.55 GPA u lt = 0.22 

L * 

a * -0.07xl0" 6 °C~ 1 « ■ 30.1xl0' 6 °C 1 

L T 

The above parameters are defined in the principal material 
system of a layer and follow the usual notation. For the case 
of temperature gradient through-the- thickness of the plate, 
results are presented for e - 0.05 in Eq. (1.1). This value of c 
means that the plate is 5 % warmer on the bottom than on the 

top. 


A continuously varying change in temperature from an ambient 
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temperature is represented by a sequence of incremental 
temperature changes. At each increment a new solution must be 
found. The search for a new solution, at each temperature 
increment, is begun by taking as an initial guess the solution 
obtained for the preceding temperature step. A total of nine 
terms are taken in the series for uo and vo in Eqs. (2.9) and 
(2 10) to approximate the inplane displacements. Typically, a 
total of four terms are taken in the series for w to approximate 
the transverse deflection. The selection of most significant 
terms for a particular laminate and the loading is generally a 
matter of trial and error Various combinations of terms are 
tried until the predominant modes for a particular case under 
consideration are identified. Transverse deflection W at the 

O 

center of plate ( normalized with respect to laminate thickness 

H) is plotted as a function of T, the change in temperature at 

* 

the midplane of the laminate ( normalized with respect to T ) . 

For plates subjected to uniform temperature, T represents the 

temperature change at every point in the plate, not just at the 

the midplane. The buckling temperature (30.6°C) of the 

four-layer square antisymmetric angle-ply (theta=45 ) plate with 

* 

fixed simple supports has been taken as the reference value T 
for antisymmetric angle-ply laminates, while for cross-ply 
laminates T* is the buckling temperature (28.6°C) of a square 
symmetric cross-ply laminate with fixed simple supports. These 
values are very close to those for an orthotropic plate (i.e. 



21 


for large number of layers) and hence are taken as the reference 
values 

Results are presented for the following two cases : 

(i) Plate subjected to uniform temperature (Figs. 3.1 to 3.42) 

(il) Plate subjected to non-uniform temperature field (Figs. 
3.43 to 3.51) 

3.2 UNIFORM TEMPERATURE RESPONSE 

Antisymmetric anele-plv laminates 

Consider the postbuckling response of an antisymmetric 

angle-ply laminate subjected to a uniform increasing temperature 

field T. Using Eq.(A1.14) of Appendix 1, it can be shown that 

for this case the thermal force N T and the thermal moments M T , 

xy x 

M T are zero. For this type of laminates, the boundary 
y 

conditions become 

Fixed simple supports ( type 1) : 

along x * 0, a 

uo = 0, vo = 0, w =0 and 

M = B (U.y + V,x + W,x W,y ) - D W,xx - D W.yy = 0 
x 1 6 * 1 

along y = 0,b 

Uo = 0, vo = 0, w = 0 and 

M = B (u.y + V,X + w , x W,y ) - D w.xx - D w.yy = 0 
y 26 “ 

(3.1) 

Sliding simple support (type 2) : 

CENTRAL LIBRARY 

I l T„ KANPUR 

MA..UHS7 
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along x =0, a 

W = 0, M x = B i6 (u,y + V.x + M , x W , y ) - D^W.xx - D^W.yy = 0 

U ° “ °* N xv A fifi (U,y + V > x + W * x W,y ) - B W.xx - B W, yy = 0 
x y o o 1 6 26 

along y =0, b 

w - 0, M = B (u.y + V , x + W , x W.y ) - D W,xx - D W,yy = 0 
y * ° 12 22 

Uo = 0, N^= A 6g (u,y + V,x + W,x W,y ) - B 16 W,XX - B^W.yy = 0 

(3.2) 

* 

It may be noted that the displacement components given by 
Eqs.(2.9) and (2.10) satisfy the essential as well as the 
natural boundary conditions. 

The results presented in Figs 3.1-3.29 are in following order : 
Figs. 3. 1-3. 3 effect of the number of layers on the transverse 
deflection of square and rectangular plates with 
lamination angle 45° and fixed simple supports 
Figs. 3. 4-3. 6 effect of the aspect ratio on the transverse 
deflection of plates with lamination 
angle 45°and fixed simple supports 

Figs. 3. 7-3. 9 effect of the number of layers on the transverse 
deflection of square and rectangular plates with 
lamination angle 22.5° and fixed simple supports 
Figs. 3.10-3.12 effect of the aspect ratio on the transverse 
deflection of plates with lamination 
angle 22.5° and fixed simple supports 
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Figs 3.13-3 21 effect of the lamination angle on the 
postbuckllng response 

Figs 3.22-3.25 effect of the number of layers and the aspect 
ratio on the transverse deflection of plates 
with sliding simple supports 

Figs. 3.26-3.29 comparison of the effect of two types of 

boundary conditions on the postbucklmg response 

In Figs. 3.1-3 3, as the temperature increases from zero, 

the plate remains flat until the temperature reaches the 

critical value required for buckling. Beyond this temperature 

the deflection increases rapidly. It is noted that the critical 

temperature at which bifurcation (buckling) occurs increases 

with the increase in the number of layers. Moreover, the 

central deflection corresponding to a given temperature load 

decreases with Increasing N. Since the coupling stiffnesses 

and B can be shown to decrease with increasing N (see 
26 

Appendix 2), it follows that membrane-twist coupling tends to 
increase the thermal lateral deflections. The difference 
between central deflections for N=4 and N=6 is very small; 
therefore, solutions corresponding to N=4,6 are close to the 
orthotropic solution. Hence, it follows that the effect of 
coupling stiffnesses is quite pronounced in N=2 case and it 
becomes less significant from N=4 onwards. For these laminates, 
the transverse displacement consists predominantly of mode (1,1) 
over the temperature range investigated; the next most 
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significant modes (1,3), (3,1) and (3,3) contribute only 

slightly to the lateral displacement. That is why, a rising 

trend in the central deflection is observed even at higher 

temperature loads, though at a reduced rate. This phenomenon is 

very clear for plates having the aspect ratio 1.5 and 2 0. As 

the aspect ratio Increases the contribution of (1,3) mode to the 

central deflection increases. This explains the trend in Figs. 

3 4-3.6, in which the central deflection of the plate having an 

aspect ratio 1.5 is less than that of plate with an aspect ratio 

2.0 at small temperature loads. But at certain critical value 
* 

of T/T, this trend is reversed. The central deflection of a 
square plate is always less when compared with the other two 
cases. 

The effect of the number of layers on the postbuckling 
response of a plate with lamination angle 22.5° can be seen in 
Figs 3.7-3. 9. In Fig. 3.7, it is noticed that beyond T/T*sl.5 
the transverse deflection of a two-layer laminate is less than 
that of a four- and a six-layer laminate. The reason for this 
is that the mode changes from (1,1) to (1,3) mode The change 
in mode shape is shown in 3-D plots (Figs 3.52 and 3.53). It is 
seen that the effect of (1,3) mode disappears as the number of 
layers increase from two to six. The effect of (1,3) mode in 
N=*6 case can be seen only at initial temperature loads when it 
produces a kink in the curve. At higher temperature loads, only 
the (1,1) mode is dominant; the (1,3) mode contributes very 
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little to the central displacements. The variation of the 
lateral displacement at the center of the plate with different 
aspect ratio is shown in Figs. 3 10-3 12. From all these 
Figures it is observed that the plate having a higher aspect 
ratio deflects more but buckles at a lower temperature. 

The central deflection and temperature relationships for 
two lamination angles considered are plotted m Figs 3.13-3.21. 

It is observed that the central deflection of a square laminate 
with lamination angle 22.5°is less when compared with the other 
case. But as the number of layers increase, the difference in 
the central plate deflection decreases and this becomes 
insignificant for N=6 (see Figs. 3.16 and 3.19). However, the 
reverse trend is observed in the case of rectangular plates. It 
is also found that the buckling temperature of a 45°laminate is 
always higher than that of a 22.5° laminate. 

In Figs. 3.22-3.25, the response of antisymmetric angle-ply 
laminates with sliding simple supports is shown. Similar to the 
case of fixed simple supports, an increase in the number of 
layers results in a decrease in the central deflection of the 
plate. Again it is due to decrease of the coupling stiffnesses 
B and B . The postbuckling behaviour of antisymmetric 
angle-ply plates for different aspect ratio is given in Figs. 
3.24-3.25. It found that the postbuckling strength decreases 
with an increase in the aspect ratio. For all these cases, the 
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mode shape (1,1) is predominant for the entire range of 
temperature considered. 

The effect of inplane boundary conditions on the 
postbuckling response can be seen in the Figs 3.26-3.29. It is 
noted that the plate with fixed simple supports buckles at a 
lower temperature. The central deflection of square plates with 
sliding simple supports is always less than that of plate with 
fixed simple supports. But, in rectangular plates having 
sliding simple supports, the out-of -plane displacement at the 
center of plate is less when compared with fixed simple supports 
only up to certain temperature load. At higher temperature 
loads out-of-plane deflections at the center of plate are less 
for fixed simple supports. In the sliding simple supports, the 
postbuckling responses consist mainly of the deflections 
associated with the coefficient w u . In the fixed simple 
supports, the plate buckles predominantly in the (1,1) mode. 
However, the deflection associated with w i3 , begin to influence 
the postbuckling response more and more at higher temperatures 
in the rectangular plates. 

Symmetric cross-ply laminates 

Using Eq. (A1.14) of Appendix 1, it can be shown that for 

this case the thermal force N^and the thermal moments M x> M y , 

M t vanish. For this type of laminate boundary conditions 
xy 


become 



Fixed simple supports: 

along x = 0, a 

uo = 0, vo = 0, w =0 and 
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M = - D W, xx - D W,yy = 0 
x 11 12 

along y = 0,b 

uo = 0 , vo = 0 , w = 0 and 


M y = " d 12 w - x>{ " D 22 w.yy = 0 


(3.3) 


The displacement components Eq. (2.9) satisfy the essential as 
well as the natural boundary conditions. 


Figs 3.30-3.32 show central deflection versus temperature 
relationships for symmetric cross-ply laminates It is noted 
that the postbuckllng strength of laminate decreases with the 
Increase in the aspect ratio of the plate The effect of the 
number of layers is shown in Figs. 3.31 and 3.32. 


Antisymmetric cross-plv laminates 

T T 

For this type of laminates N and M are zero. Therefore, 

xy xy 

the boundary conditions become 


Fixed simple supports: 

along x ** 0, a 

uo = 0, vo * 0, w =0 and 
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V B ,, (u - X + l/2(w,x) 2 ) - D n W,xx - D^W.yy - M*= 0 

along y = 0,b 

uo*0, vo = 0, w = 0 and 

M = 8 (v.y * »/2(w,y) 2 ) - D W.xx - D W,yy - M T = 0 

y 12 22 y 

(3.4) 


Sliding simple supports: 
along x * 0, a 

w * 0, M * B (u,x + i/2(w, x ) 2 ) - D w i xx - D W,yy - M T = 0 

X 1 1 11 12 x 

Uo = 0. N = A (u, y + V,x + W , x W,y ) = 0 

x y 6 6 

along y * 0, b 

2 x 

W = 0* M « B (v»y + 1/2 (w,y) ) “ D W,xx - D W,yy-M = 0 

y 22 12 22 y 

Vo *= 0* N = A (u,y + V,x + W, x W t y ) - 0 

xy 66 

(3.5) 

In this case the assumed displacement components, expressed by 
Eqs. (2.9) and (2.10), satisfy the essential boundary conditions 
but not the natural boundary conditions. However, as in 
previous cases, nine terms in u and v displacement series and 
four terms in the w series are considered. 

Figs 3.33-3.36 give transverse central deflections for 
square as well as rectangular cross-ply laminates. Unlike the 
cases considered earlier in which the simply supported laminates 
remain flat until the bifurcation (buckling) occurs, the 
antisymmetric cross-ply laminate deflects laterally immediately 
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upon loading. This unique behaviour of cross-ply can not be 

noticed by linear analysis (Thangaratnam, et al 1988). This 

behaviour is due to the presence of coupling stiffnesses B^and 

B 22 which couple normal moments with inplane compression. Fig 

3 33 shows the effect of the number of layers on the transverse 

deflection of cross-ply laminates. It is observed that there is 

hardly any difference between the responses of four- and 

six-layer laminates Since the coupling stiffnesses B^and 

B can be shown to increase with decreasing N, it follows that 
2 2 

B and B tend to Increase deflections. In rectangular plates 
1 1 22 

(a/b = 2.0, see Fig. 3.36), a two-layer laminate buckles in 
(1,1) mode but (1,3) mode begins to influence the response more 
and more as w approaches the magnitude of w^. At a specific 
temperature the central deflection of a two-layer laminate 
becomes less than that of a four-layer laminate. Here at this 
point it can be noticed that in the of four-layer laminate, the 
curve has a rising trend due to the presence of mode (1,1) only. 
Figs 3.34 and 3.35 show the effect of the aspect ratio on the 
central deflection of the plate. 

Figs. 3.37-3.40 depict the nonlinear response of cross-ply 
laminates with sliding simple supports (type 2 boundary 
condition). Comparative study of the influence of two different 
boundary conditions is given in Figs. 3.41 and 3.42. It is 

with the fixed simple supports deflect 


observed that plates 
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out-of-plane more than plates with the sliding simple supports. 


Angle-ply Vs . Cross-plv laminates 

On the basis of the study carried out for antisymmetric 
angle-ply and cross-ply laminates the following points are 
observed : 

(i) The most interesting difference in the behaviour of two 
types of laminates is that the cross-ply laminates deflect 
out-of-plane as soon as thermal loads are applied, whereas 
angle-ply laminates remain flat till a critical temperature is 
reached after which the deflection increases rapidly with rise 
in temperature. In other words, bifurcation type of instability 
(buckling) is not observed in the case of antisymmetric 
cross-ply laminates. 

£il) Out-of-plane deflection at the center of plate does not 
depend significantly on the number of lay-ups beyond N=4 in both 
angle- and cross-ply laminates (Figs. 3.1 and 3.33). 

(iii) The central deflection increases with the increase in the 
aspect ratio in both the cases considered but mode shape may 
influence the central deflection of the plate (see Figs. 3.5 and 
3.39). 

(iv) The effect of inplane boundary conditions on the transverse 
deflection is much more pronounced in the case of cross-ply 

be observed that that the plates with 


laminates, where it can 
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the sliding simple supports are stiffer compared to the plates 
with the fixed simple supports (Figs 3.41 and 3.42). 
Ant isymmetric angle-ply laminates do not show any significant 
difference in the lateral deflections for the two boundary 
conditions considered ( Figs. 3.26 and 3.29). 

3.3 NON-UNIFORM TEMPERATURE RESPONSE 

Figs 3.43 to 3.51 show the response of antisymmetric 
angle-ply laminates in the presence of thermal gradient 
through-the-thickness of the plate. It is found that the plate 
undergoes lateral deflection as soon as the temperature is 
increased above the ambient temperature. However, for both sets 
of boundary conditions the deflections prior to the critical 
buckling temperature are very small and at only slightly higher 
temperature, the response becomes very close to the 
corresponding postbuckllng response in the uniform temperature 
case. The effect of the number of layers on the out-of -plane 
deflections in the presence of slight thermal gradient is also 
studied. It is noticed that the response with higher number of 
layers is even closer to the postbuckling response of the 
earlier case. Therefore, it appears reasonable to conclude that 
a slight temperature gradient through-the-thickness of plate 
does not make much difference in the nonlinear response of the 
plate. A similar observation has also been reported recently by 
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Meyers and Heyer (1992). 

It is further observed that the response of the plate with 
different fiber orientations in the presence of the temperature 
gradient is same as when the temperature is uniform, except that 
initial out-of-plane deflections are observed in the 
temperature-gradient case (Figs. 3 43-3.48) Again, the central 
deflection of the square plate with lamination angle 22.5° is 
less when compared with the case of lamination angle 45° But 
as the number of lay-ups N increase, the difference in the 
out-of-plane displacement at the center of the plate decreases 
for the two cases and this difference become insignificant for 
N=6 (Figs. 3.45 and 3.48). The effect of inplane boundary 
conditions can be seen in Figs. 3.43-3.45 and Figs. 3.49-3.51. 
Here it Is again noted that the central deflection of a square 
plate with sliding simple supports (type 2) is always less than 
that of a plate with fixed simple supports (type 1). Although 
the results for rectangular plates are not given here but it has 
been verified that their general response is similar to that for 
uniform temperature case except for the out-of-plane deflections 
setting in as soon as thermal loads are applied. 
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Fig. *• *3 Contor d#f Wet ion of anlisymmotric 
angW-pty »ubj*ct*d l© uniform Wmporatur* 
Boundary condition*: Typo 1 



T/r 


Fig* B« 14 C#ni*r dofWction of ant vsyrrmit trie 
arvgW-ply subjoctsd to uniform i*mp*raiur* 
Boundary conditions: Typo 1 



Boundary condUtons: Typ# 1 




38 



Fig* a » Gontor dofloction of antvoymmetric 
ar»g W~ply *ub joclod to uniform tomporaturo 
Boundary condition#: Typ* * 



Fig. 3. 17 Gontor dofloction of antisymmetric 
anglo-ply subjected to uniform temperature 
Boundary condition*: Typo 1 



Fig* 3. 18 Center deflection of antisymmetric 
angle-ply •ubjected to uniform temperature 
Boundary condition*: Typo l 
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CONCLUSIONS 


4 1 CONCLUDING REMARKS • 

The method of incremental thermal loads used in this 
investigation provides a means for studying moderately large 
deformation response of angle- and cross-ply laminated plates 
Two types of inplane boundary conditions have been considered. 
The effect of imperfection in temperature field in the form of a 
slight temperature gradient through- the-thickness of the plate, 
has been considered. 

On the basis of the limited study carried out the following 
conclusions can be drawn: 

(i) The thermal postouckling characteristics are 
significantly influenced by the number of layers. The effect of 
coupling stiffnesses B^and B^ is quite pronounced in a two-layer 
case and It becomes less significant from the four-layer case 
onwards. The change in mode shape also plays very important role 
in determining the lateral displacement of the plate. 

(li)The thermal postbuckling strength decreases with the 
increase in the plate aspect ratio for angle-ply laminates and 
cross-ply laminated plates, 

(ill) The magnitude of the transverse deflection at the 
center of the plate gets substantially changed by changing the 


lamination angle 0. 



(ivl Support conditions can change the deflection response of the 
laminates In general, plates with fixed simple supports deflect 
more than plates with sliding simple supports. 

(v) Antisymmetric cross-ply laminates deflect out-of-plane 

immediately upon the application of thermal loads, i e the 

classical buckling is precluded in this case. This unique 

behaviour is due to presence of B and B terms in the coupling 

1 1 22 

nM t r i x 

(vi ) The presence of small imperfection in the temperature, m the 
for m of through-the-thickness gradient does not greatly affect 
nonlinear behaviour of laminates. 

4.2 RECOMMENDATIONS FOR FUTURE WORK * 

The formulation developed does not include transverse shear 
deformation effects. As a/H ratio becomes less than or equal to 
20, the* shear deformation effects become increasingly important. 
This aspect can be included in the analysis Most of the studies 
carried out recently consider only lateral defections of the plate 
even at very high temperature loads. At that stage excessive 
stresses may cause delamination, cracks etc Therefore, it is 
equally important to include stress analysis in the investigation. 
The combined effect of mechanic and thermal loads on the 
postbuck ling response of laminate is also an important topic for 


study. 
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APPENDIX 1 


1.1 STRESS-STRAIN RELATIONS FOR A LAMINA . 

If the material axes are aligned with the natural body axes 
for the problem, then such lamina is called specially orthotropic 
lamina. Stress strain relations for such lamina under plane 
stress can be written as (Agarwal 8. Broutman 1980) 
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where L and T subscripts refer to the longitudinal direction (L) 
and the transverse direction (T), respectively. [Q1 Is the 
reduced stiffness matrix and can be expressed in terms of 
engineering constants as under 
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(A1.2) 


Here E E are the elastic moduli in L and T directions (Fig 
L * T 

Al.l), respectively, G is the shear modulus in the L-T plane, v 

J LT 

is the Poisson’ s ratio for transverse strain in the T-direction 
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when ».tti“sscd in the L-direction 


d 



T 


Fig. Al.l Orthotropic lamina oriented along its 
principal coordinate axes. 

The pi incipal directions of orthotropy often do not coincide 
with coordinate directions that are geometrically natural to the 
solution of the problem (Fig. A1.2). Stress strain relationship 
for such lamina under plane stress are 



where the transformed reduced stiffness Qij can be written In 
terms of reduced stiffness Qij in the relieving manner 
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where C = Cos9 , S = Sine and 0 is the angle from X-axis to 
principle orthotropic axis (L-axis) as shown in Fig.(A1.2) 
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Fig. A 1 . 2 Ortho tropic lamina with its principal 
axes oriented at angle 0 with reference coordinate 
a x e s . 


Thus, in body coordinates, even an orthotropic lamina appears 
to be anisotropic. However, because such a lamina does have 
orthotropic characteristics in principal material directions, it 
is called generally orthotropic because it can be represented by 
stress-strain relations given by Eq. (A1.3). The Eq.(A1.3) can be 
thought of as stress-strain relations for the k th layer of a 
multilayered laminate. Thus it can be written as 
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1.2 STRAIN VARIATION IN A LAMINATE : 

The laminate is presumed to consist of perfectly bonded 
laminae. The bonds are infinitesimally small as well as 
non-shear -deformable; that is, the displacements are continuous 
across lamina boundary so that no lamina can slip relative to 
another. The Kirchhoff’s hypothesis implies a linear variation of 
strain through the laminate thickness. The strain in the k th 
layer can be expressed as 
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where the middle surface strains are 
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and the middle surface curvature can be expressed as 
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Here uo and vo are the middle 


surface displacements in the x- and 
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y-dirt?< t ions, respectively, and w is the lateral displacement. 

1 3 RESULTANT FORCES AND MOMENTS OF A LAMINATE. 

The resultant forces and moments acting on the cross-section 
of the laminate are defined as follows 
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These force (N , N , N ) and moment (M , M , M ) resultants 
x y xy x y xy 

for a laminate can be related to the middle surface displacement 

j j* j 

components (uo, vo, w), thermal forces (N , N , N ) and thermal 

x y *y 

moments (M*, M*. by using expression (Al.ll) and (A1.12) 
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Here A (j , and D (j are called the extensional, the coupling and 
the bending stiffnesses, respectively, and are defined in terms of 
*' or tlle incIiv *dual layers k * 1, 2 N : 

^ r ^ 

(A U' V D ij ) 3 Y* C1, Z ’ z2)( ^JJ } lc d2 (A1.13) 

k=l z 

k-1 

(i, j) = (1, 2, 6) 


The limits z k and z^are shown in Fig. A1.3. The thermal forces 
and moments appearing in Eq.(Al.ll) and Eq.(A1.12) are defined as 



(A1.14) 


where a , a are the coefficients of thermal expansion in X- 
x y 

and Y-directlon respectively, and a is an apparent coefficient 
of thermal shear. 



Fig. /\1«3 Geometry of an N-tayred laminate (Jones, 
1975) 


APPENDIX 2 


1. Special cases of laminate stiffnesses 


1 1 Antisymmetric Angle-ply Laminate 
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1.2 Symmetric Angle-ply Laminate 
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1.5 Antisymmetric Cross-ply laminate 
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It can be seen from Eqs. (A2. 1 )-(A2. 4) that as the number of 

layers increases (i) the coupling stiffnesses tend to zero (ii) 

the extensional stiffnesses A and A , in the case of symmetric 

16 26 

angle-ply laminate, tend to zero 

Hence, the laminate behaves as a single specially orthotropic 
layer since B , A , A ,D and D are zero. Hence the solution 

Ij 16 26 16 26 

corresponding to N = « is termed as the orthotropic solution. 



